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The freeze out of the expanding systems, created in relativistic heavy ion collisions, 
will be discussed. We combine kinetic freeze out equations with Bjorken type system 
expansion into a unified model. Such a model is a more physical generalization of the 
earlier simplified non-expanding freeze out models. We shall see that the basic freeze out 
features, pointed out in the earlier works, are not smeared out by the expansion. 



1. Model for the simultaneously expanding and freezing out 
fireball 

In this paper we present a simple freeze out (FO) model, which describes the FO 
of particles from a Bjorken expanding fireball ^. The important feature of the 
proposed scenario is that physical freeze out is completely finished in a finite time, 
which can be varied from (FO hypersurface) to oo. In the other words our FO 
happens in a layer, i.e. in a domain restricted by two parallel hypersurfaces r = ri 
and r = Ti + L, where r is the proper time variable and L is the maximal duration 
of the FO process. 

The present model simultaneously describes the gradual FO and the expansion 
of the system, and thus it is a more realistic ext ension of the oversimplified FO 
models, which did not include system expansion I 2 | 3 | 4 | 5 | Ref. ^ authors have 
also adopted kinetic gradual FO model to Bjorken geometry, but combined it with 
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Bjorken expansion on the consequent, not on the paraUel basis: system expands 
according to Bjorken hydro scenario, but when it reaches beginning of the FO 
process system stops expansion and gradually freezes out in a fixed volume. It was 
shown in 1^ that although such a model allows to obtain analytical results, it is not 
physical, the simultaneous modeling of expansion and freeze out is required in order 
to avoid decreasing of the total entropy. And now we propose such a generalized 
model 1^. 

We start with introducing two components of the distribution function, /: the 
interacting, and the frozen out, ones, (/ = f^ + f^), then, correspondingly we 
will have two components of the energy density and baryon density. Details of the 
derivation can be found in Ref. ^ and here we only give the final coupled system of 
equations for interacting and free components describing the change in the particle 
density and energy density: 

dr T Tpo L + Ti — T ' dr 

de^ e* L dn^ 

dr T Tpo L + Ti — r ' dr 

All together we have the following simple model describing the evolution of the 
fireball created in relativistic heavy ion collision. 
Initial state: r = Tq , 

e(T) = eo , n(r) = no . (3) 
Phase I, Pure Bjorken hydrodynamics: to < r < ti , 




where P — CqC is the equation of state (EOS). 

Phase II, Bjorken expansion and gradual FO: t\<t<t\+L , 

e\r)=eo[-) [^^) ' " = "° (v) j '^'^ 

The difference with respect to the pure Bjorken solution eqs. ^ is the multiplier 
describing the gradual FO of the system. We also see that the interacting component 
will vanish when we reach the end of the FO layer, i.e. r — > L + ti. 

Knowing e*(T) and equation of state (EoS) we can calculate the temperature of 
the interacting component, Ii(r), as a function of the proper time. Due to symme- 
try of the system, uf (t) — u^{tq) — (1,0, 0, 0). Thus, we have complete knowledge 
about the evolution of the interacting component, /'(t), which is a thermal distri- 
bution with given Ti(T), ri*(T), u^{t). 

However, what we have to calculate is the free component, which is the source 
of the observables. Eqs. ([2]) give us the evolution of the e/ and n/, and one can 
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= 5.0 fm/c, XpQ=0.5 fm/c 




X, [fm] 

Fig. 1. Evolution of temperature of the interacting matter for different FO layers. Ti{To = 
0.05 fm) = 820 MeV, Tpo = 180 MeV. "No Freeze Out" means that we used standard Bjorken 
hydrodynamics even in phase II. 



easily check that these two equations are equivalent with the following equation for 
the free component: 



dr T tfo \L + Ti — T 

The measured post FO spectra are given by f-f{L + Ti). 

2. Results from the model 

Aiming for a qualitative illustration of the FO process we show below the results 
for the ideal massive pion gas with Jiittner equilibrated distribution 



where the degeneracy of pion is .g = 3, while the baryon chemical potential in case 
of pions is zero. 

Contrary to the illustrative example in here we do not neglect the pion mass. 
During FO the temperature of the interacting component decreases to zero, so at 
late stages of the FO process this new calculation is better justified. We will see 
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below that Ti falls below TO^r quite soon, and so the Jiittner distribution is a good 
approximation of the proper Bose pion distribution. 
For our system we have the following EoS: 

where Kn is Bessel function of the second kind, and a — m/Ti. 

The first eq. of system ([T]), gives the following equation for the evolution of the 
temperature of the interacting component: 

<m ^ T, 4TfK2{a) + mT,Ki{a) 

dr T 12T^K2(a) + 5mTiKi{a) + m^Ko{a) 



— 1 

tfo \L 



?>TfK2{a)+mT,Ki{a) 



(9) 



Ti-T J l2TfK2{a) + 5mTiKi{a) + m^Ko^a) ' 

Furthermore, we have used the following values of the parameters: ijr = 4.38, A^y — 
TrR\^, where Rau = 7.685 fm is the Au radius, tq = 0.05 /m, T,(to) = 820 MeV, 
Ti = 5 fm, what leads to Ti(Ti) = Tpo — 180 MeV, and Tpo = 0.5 fm. During 
the pure Bjorken case the evolution of the temperature is govern by a simplified eq. 
(I9|), without the second freeze out term on r.h.s. 

In Fig. [l]we present the evolution of the temperature of the interacting matter, 
Ti{T), for different values of FO time L. 
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Fig. 2. Final post FO distribution for different FO layers as a function of the momentum in the 
FO direction, p = in our case (p" = p^ = 0). The initial conditions are specified in the text. 



As it was already shown in Ref.^ ^, the final post FO particle distributions, 
shown on Fig. [21 are non-equilibrated distributions, which deviate from thermal 
ones particularly in the low momentum region. By introducing and varying the 
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thickness of the FO layer, L, we are strongly affecting the evolution of the interacting 
component, see Fig.[Tl but the final post FO distribution shows strong universality: 
for the FO layers with a thickness of several Tpo post FO distribution already 
looks very close to that for an infinitely long FO calculations, see Fig. [2] left plot. 
Differences can be observed only for the very small momenta, as shown in Fig. [2] 
right plot. So, the inclusion of the expansion into our consideration does not smear 
out this very important feature of the gradual FO. 

It is important to always check the non-decreasing entropy condition EH to 
see whether such a process is physically possible. Figs. [3] present the evolution of 
the total entropy, S{t), calculated based on the full distribution function, /(p) = 
r(p) + /^(p): 



1 - In 
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fir) 



S{t) = s{t)V{t) . 



(10) 



During pure Bjorken phase total entropy remains constant, as expected, but during 
phase II it constantly increases until FO is finished. 
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Fig. 3. Evolution of the total entropy for different FO layers. The initial conditions are specified 
in the text. 



From this figure we can make an important conclusion, that gradual freeze 
out produces entropy. For the late stages of the FO this can be approximated 
analytically. 



ds{T) 
dT 



,3 df{T) 



dT 
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dfjr) 
dr 



(11) 



At the late stages of the FO contribution of the interacting component can be 
neglected, ao f ~ f-^ , and, from eq. ([6|), ^ « — So for the late stages of the 
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reaction we have the following equation for the entropy evolution: 

ds{T) s{t) (t) 
dr T T 

And since in the Bjorken geometry V{t) ~ r we obtain: 

dS_d{sV)^ sV ni'V sV _ Nf 

dr dr t t t t ' 

where Nf is total number of frozen out particles. Thus, we see that total entropy 
increases during the simultaneous expansion and gradual FO, and at the very end 
of the FO process it increases logarithmically. 



(12) 



(13) 



3. Conclusions 

In this paper we presented a gradual FO model including a Bjorken like expansion, 
being an extension to the older versions 2 3 4 5^ which allowed us to study FO in a 
layer of any thickness, L, from to oo. Another important feature of the proposed 
model is that it connects the pre FO hydrodynamical quantities, like energy density, 
e, baryon density, n, with post FO distribution function in a relatively simple way, 
and furthermore allows analytical analysis for simple systems, like massless pion 

The results show that the inclusion of the expansion into FO model, although 
strongly affects the evolution of the interacting component, does not sm ear ou t the 
universality of the final post FO distribution, observed already in Refs. I3HE1. 

Another important conclusion of this work, stressing once again the importance 
to always check the non-decreasing entropy condition IHIII]^ jg that long gradual 
freeze out may produce substantial amount of entropy, as shown on Fig. O 
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